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We report on experiments with superconducting metamaterials containing Josephson junctions. 
In these structures, split-ring resonators used in conventional metamaterials are replaced by su- 
perconducting loops that are interrupted by Josephson junctions, so called rf-SQUIDs. Like the 
split-ring resonators, these elements can be seen as LC-resonators that couple to the magnetic field. 
The advantage of superconducting thin- film metamaterials is that, due to the tunable intrinsic 
inductance of the Josephson junction, the resonance frequency of the rf-SQUID can be changed 
by applying an external dc magnetic field. We present experimental results that demonstrate the 
tunability of the resonance frequency of these devices. 



One of the main limitations of many metamaterial de- 
signs is the restriction of their usability to a narrow fre- 
quency band. By using resonant elements in order to 
achieve the desired effects (such as a negative magnetic 
permeability fi r ) the meta- atoms are usually made for a 
fixed operation frequency range by design. To circumvent 
this limitation, several approaches have been developed 
[l| , many of which are based on the idea of introducing a 
nonlinear element into the resonant circuit, thus making 
it tunable [H, Q . Superconductors in general and Joseph- 
son junctions in particular are ideal constituents of tun- 
able meta-atoms due to the strong dependence of their 
inductance on external parameters such as temperature 
and magnetic field. The tunability by varying tempera- 
ture has already been demonstrated with superconduct- 
ing split-ring resonators (SRRs)0, |]|. In this work, we 
move one step further and experimentally demonstrate 
tunability via magnetic field by replacing the SRR by a 
superconducting quantum interference device (SQUID). 
This approach has previously been suggested in the the- 
oretical work by Lazarides and Tsironis J6j and further 
investigated by Maimistov and Gabitov [7j. 

The key element used in this work is the Josephson 
junction which consists of two superconductors separated 
by a thin insulating layer. The current and voltage over 
such a J unction is described by the two Josephson rela- 
tions |8l-[lQ|. The Josephson junction with a critical cur- 
rent I c can be treated as a nonlinear network element. A 
simple description is provided by a small signal approx- 
imation. A small variation Sep of the phase difference 
if between the superconducting wave functions on both 



sides of the junction around a constant value tp' is trans- 
lated into a small variation of the current I = I c sin (p'+SI 
flowing through it. In this way, current and voltage are 
related to each other by [H, [Tof : 

27r/ c cos(^ / dt JVV ^ ' dt w 

This equation takes the form characteristic to an in- 
ductive element. That is, the voltage over the element is 
proportional to the time derivative of the current through 
it. In this model, the value of the so-called Josephson in- 
ductance Lj depends on the constant part of the phase 
difference across the junction and therefore on the bias 
current through it. Thus, it works as a tunable inductor. 

A circuit consisting of a Josephson junction placed into 
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FIG. 1: (a) A superconducting loop interrupted by a Joseph- 
son junction can be decomposed into a (b) network model 
using the restively and capacitively shunted junction model 
(RCSJ) pjj. The red cross indicates the Josephson junction 
for which voltage and current relation are given by the Joseph- 
son relations, (c) Small signal network model of the SQUID. 
Lgeo presents the geometric inductance of the superconduct- 
ing loop. Lj is the Josephson inductance, see Eq. (pQ). R and 
C are the normal resistance and capacitance of the Josephson 
junction, respectively. In our sample, the junction is shunted 
by a large additional capacitance to decrease the resonance 
frequency, which is added in parallel to C. 
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a superconducting loop (cf. Fig.QJi) is a single junc- 
tion (rf-) SQUID. In this case, magnetic flux quantization 
leads to the relation 

(f^ + L geo I = $ ext . (2) 

The electric circuit equivalent of an rf- SQUID is de- 
picted in Fig.[T|3. In addition to the Josephson relations, 
one has to account for the loop inductance, the capac- 
itance between the superconducting electrodes and the 
(resistive) single electron tunneling current. 

Once an external magnetic flux 3> e xt(£) = ^eo + 
<£ ea cos(ut) is applied to the system and the driving ac 
amplitude <£ ea <C I c /L geo is small, the small signal ap- 
proximation is valid. The total inductance L tot comprises 
the geometric inductance L geo and the tunable Josephson 
inductance Lj(<J> e o), see Fig. [It. The resonance frequency 
of the system is given as a function of the external bias 
flux $ e0 

M$e0) = (Ltot(*e )C)-i • (3) 

In the present experiment, the additional kinetic induc- 
tance [9] of the superconducting ring is negligible com- 
pared to its geometric inductance. It should be noted, 
though, that close to the critical temperature of the su- 
perconductor, the kinetic inductance will increase signif- 
icantly and influence the resonance frequency. Due to 
this fact and the temperature dependence of the criti- 
cal current of the junction, the rf-SQUID also shows an 
intrinsic tunability by temperature which, however, will 
not be detailed here. 

The approximation of the Josephson junction by a tun- 
able inductor works well within the mentioned small sig- 
nal limit. However, in order to obtain a complete descrip- 
tion, the Josephson nonlinearity and loss terms have to be 
included. By applying Kirchhoff 's laws, the two Joseph- 
son relations and Eq. (|2]) to the circuit depicted in Fig.[T]3, 
a second order nonlinear differential equation is derived: 




FIG. 2: Sketch of a single rf-SQUID embedded into a copla- 
nar transmission line. The darker areas represent the Nb elec- 
trodes. The inset on the right shows an optical microscope 
image of a single rf-SQUID. 



The time average of the solution of this equation yields 
the complex ac flux susceptibility of the single SQUID 
which is x</> = ($/$ext) t — 1- There is a strong fre- 
quency dependence of real and imaginary components 
of Xcf> around the resonance frequency ujo. 

For an effective medium consisting of these elements, 
the magnetic susceptibility is given by Xmag = ^X</>- The 
coefficient F depends on the arrangement of the meta- 
atoms in the medium. Provided that the coupling be- 
tween the SQUIDs is small (i.e. F <C 1) [y], the relative 
permeability of the medium is 

Mr = 1 + Xmag = 1 + ^X</>- (5) 

In this work, we experimentally verify the microwave 
properties of a prototype circuit representing a single su- 
perconducting meta-atom of such a medium. We have 
designed and fabricated a circuit consisting of a single 
rf-SQUID placed inside the gap of a coplanar waveguide 
(CPW). In this setting, the SQUID loop is perpendicular 
to the magnetic field, as can be seen in Fig. [2^, and thus 
couples to the quasi-TEM mode of the CPW. The arrows 
denoted by E, H and S indicate the orientation of the 
electric and magnetic field components and the Poynting 
vector of the incoming wave, respectively. An external 
magnetic field is applied perpendicular to the substrate 
plane either through an external coil or by driving a dc 
current through the central strip of the transmission line. 
In either case, the SQUID loop is biased with a dc mag- 
netic flux ^> e o thus making the resonance frequency tun- 
able in situ according to Eq. ((3]). 

The SQUID loop used in the presented measurement 
is made from niobium (which has a superconducting 
transition temperature of T c = 9.2 K) and includes a 
Nb/A10 x /Nb circular Josephson tunnel junction with a 
radius of 0.5nm and a critical current I c ~ 4.5 \iA . The 
geometric loop inductance L geo = 63.5 pH is comparable 
to the zero-flux Josephson inductance Lj (0) = 74 pH. By 
shunting the Josephson junction with an additional large 
capacitor C s h U nt = 1.7 pF we designed the resonance fre- 
quency to be tunable between about 21 GHz and 5 GHz, 
for an infmitesimalfy small <fr e8i . The measurements were 
performed using a cryogenic amplifier and a vector net- 
work analyzer while the sample was immersed in liquid 
4 He at a temperature of 4.2 K. 

The experimental results for low-power transmission 
presented in Fig.[3fci demonstrate the frequency tunability 
of the resonance dip. The decreased transmission is illus- 
trated by the darker color in the plot. The resonance fre- 
quency changes periodically as a function of the applied 
magnetic flux <I> e o with a period of ^>o- The dashed blue 
line represents a fit of the small signal model of Eq. ((3]). 
As discussed earlier, the model should yield good fit re- 
sults only for very small values of the driving amplitude. 
The power used in this measurement was -85 dBm which 
corresponds to an ac flux amplitude of <£ ea w <I>o/1600 
and can be safely assumed to be small enough for the 
approximation to hold. The values, obtained from the fit 
for Lg eo and C agree within 5% with the values expected 
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FIG. 3: (a) Measured microwave transmission at a power 
of P ~ —85 dBm at the sample, that corresponds to <J> e a ~ 
<3>o/1600. The darker colors correspond to a reduced transmis- 
sion through the SQUID-loaded waveguide. The blue dashed 
line represents the analytical solution extracted from the small 
signal model (Fig.[T]:). It is shown for negative flux values only 
in order not to obscure the data, (b) Transmission through a 
single SQUID simulated numerically by solving Eq. (4|). 



pie, analytical description, we approximated the system 
by treating it as a homogeneous, cylindrically symmetric 
problem, similar to a coaxial cable. The unit cell of the 
waveguide, of length /, contains one SQUID which de- 
fines the interaction with the magnetic field component. 
Under these assumptions 
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Here, A sq is the area of the SQUID, r2 and T\ are the 
outer and inner radius of the gap, as shown in Fig. [2^. 

Using Eqs. ([5]) and (|6]), we obtain the effective perme- 
ability fi r and use it to modify the inductance and thus 
the impedance of the CPW section. This, in turn, allows 
us to calculate the scattering matrix coefficient S21. The 
result of this calculation, depicted in Fig. [3b, is in excel- 
lent agreement with the measured data. As expected, 
at so small applied power level, the numerical treatment 
adds very little additional information to the analytical 
solution of Eq. ([3]) shown in Fig. [3^ by the dashed blue 
line. However, at higher power levels, the small signal 
approximation of Eq. (pQ) is no longer valid and the nu- 
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for the fabrication procedure while the value of the criti- 
cal current I c is about 6% less than the one measured on 
test junctions on the same chip. 

In order to estimate the magnetic susceptibility and 
thereby the relative permeability, Eq. (j4]) was solved nu- 
merically. The parameters L geo and C were taken from 
the fit mentioned in the paragraph above. The value for 
the shunting resistance R was taken from the current- 
voltage characteristics measured for a test junction on 
the same chip. Additionally, the loss tangent of the 
Nb2C>5 dielectric of the shunt capacitor was assumed to 
be tan 5 = 1.0- 10 -3 pjj. To calculate the scattering char- 
acteristics versus dc magnetic flux <l> e o, we employed the 
effective medium approach by calculating the effective 
fjL r of a SQUID-loaded section of the CPW. Analogous 
to the procedure in Ref. [6|, our task was to find a suit- 
able value for F. Since the actual field distribution in a 
superconducting coplanar waveguide does not have a sim- 
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FIG. 4: Relative (a) magnitude and (b) phase of the trans- 
mission (S21) through the SQUID-loaded coplanar waveguide 
at a bias of <3> e o = 0.77 $o- The measured data (power 
P « —110 dBm) are shown in black, the simulated curve 
($ ea = $ /2 15 ) in blue. 



merical treatment will be relevant. 

The frequency-dependent properties of a metamaterial 
can be characterized by relative magnitude and phase of 
the complex transmission coefficient S^i- Measurement 
and numerical calculation of the transmission magnitude 
and phase through the SQUID-loaded coplanar waveg- 
uide (cf. Fig. [4^1 and HJd) are in good qualitative agree- 
ment. Detailed quantitative comparison in our case is 
hindered by difficulty of making precise calibration mea- 
surements at cryogenic temperatures. 

In conclusion, we have experimentally and theoreti- 
cally investigated a prototype of a magnetically coupled, 
tunable meta-atom made of a thin-film superconducting 
circuit. The resonance frequency of a superconducting 
loop containing a Josephson junction is tunable over a 
wide frequency range both by magnetic field and tem- 
perature. Therefore, the SQUID metamaterial appears 
as an appealing alternative to the widespread normal- 



metal structures employing split-ring resonators. Further 
studies of SQUID properties at higher powers, i.e. inves- 
tigating its nonlinear properties, as well as experiments 
on ID and 2D-arrays of such elements are presently being 
conducted. 
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